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Abstract 

A model describing the Anderson impurity in the Bardeen-Cooper-Schriffer 
superconductor is proven to exhibit hidden integrability and is diagonalized 
exactly by the Bethe ansatz. 
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The basic theoretical models describing magnetic impurities in nonmagnetic normal met- 
als, such as the s-d (Kondo) model, the Anderson model, etc., are integrable under two 
additional conditions |1|]: (i) an electron-impurity coupling is assumed to be energy inde- 
pendent, and (ii) a band electron dispersion s(k) can be linearized around the Fermi level, 
e(k) ~ vp(k — kp), where kp and Vp are the Fermi momentum and velocity, respectively 
Only under these conditions, both an electron-impurity scattering and an effective electron- 
electron coupling are described in terms of discontinuous jumps in the Bethe ansatz wave 
functions. Therefore a linear dispersion of particles and a pointlike particle-impurity cou- 
pling are considered now as the necessary mathematical conditions for integrability of the 
"impurity" models. Because a carrier dispersion in superconducting metals cannot be lin- 
earized, the "linear dispersion approximation" (LDA) is clear to eliminate superconductivity 
from an exact analysis of the behaviour of magnetic alloys Q. 

It has recently been discerned that LDA is not necessary for an exact diagonalization 
of the basic impurity models of quantum optics, describing a system of Bose particles with 
a nonlinear dispersion coupled to two-level atoms ||. Making use of some mathematical 
analogies between "magnetic" and "optical" models, it can be shown || that the degenerate 
Anderson model with a nonlinear band electron dispersion also exhibits hidden integrability 
H and is diagonalized exactly by the Bethe ansatz. One of the most exciting applications of 
the approach developed could be an exact treatment of the superconductivity problem 
in the presence of magnetic impurities. Since the pioneering work of Abrikosov and Gor'kov 
H, the problem has been the subject of many theoretical and experimental studies || but 
still remains one of the most intriguing puzzles of modern physics. Therefore an extension 
of the Bethe ansatz technique to superconducting magnetic models accounting for a gap 
dispersion of charge carriers is of fundamental physical interest. 

In the present letter, we report first results for a model describing the Anderson im- 
purity placed within the Bardeen-Cooper-Schriffer (BCS) superconductor. To diagonalize 
the model Hamiltonian, we introduce auxiliary particles and show that the multiparticle 
scattering process is factorized into two-particle ones, despite a nonlocal effective particle- 
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impurity coupling. The continuous multielectron wave functions result from an integral 
"dressing" of the discontinuous Bethe ansatz wave functions of auxiliary particles, the in- 
formation about the electron dispersion being contained into a dressing function. Imposing 
the periodic boundary conditions on the multielectron wave functions, we derive the Bethe 
ansatz equations (BAE) for the cases of a rare-earth and a transition metal impurity with 
infinitely large Coulomb repulsion on the impurity orbital. 

To derive the model under consideration, one can start with the Hamiltonian involving 
the Hamiltonians of the BCS and Anderson models, 

<r=U k k 

+ J2J2 v k ( a ka d <r + 4 a ka) + e d ^2 d x a d a + Ud\d^d\di. (1) 

cr=U k V ' a=U 

The Fermi operator a ko . creates a conduction band electron with the momentum k, spin 
a =|, | and the energy e k = e k — ep, where e k and ep are the kinetic and Fermi energies. 
An electron localized on the impurity <i-orbital with energy is described by the Fermi 
operators d a . The third term of Eq. (1) represents the hybridization of a band and a d- 
level electrons with the matrix element t> k , while the Coulomb repulsion on the impurity 
orbital is described by the last term. Two first terms of Eq. (1) are the standard BCS 
model, where the gap A is assumed to result from the Cooper pairing phenomenon. For our 
purposes, it is convenient to treat A and ep as some free parameters, which have to be found 
self-consistently, and to start our analysis with the bare vacuum of the model defined by 
o ko -|0) = d a \0) = 0. In the normal state (A = 0), the model (1) is reduced to the Anderson 
model, which has been diagonalized by Wiegmann || within the LDA. 

The BCS part of the model Hamiltonian is diagonalized by the Bogoliubov-Valatin uni- 
tary transform j7| to give 

#bcs = Y.J2 ^kC kCT c k(T , (2a) 

a- k 

where co> k = —J~£^ + \ A| 2 , for k = |k| < kp and u; k = \j£l L + | A| 2 for k > kp. In terms of 
new Fermi particles (which we will also call "electrons"), the hybridization term takes the 
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form 

y = *££(4ok, + 4,A)> (2b) 

a k 

where we have set k = kp both in the hybridization matrix element, v = v(kp), and in the 
coefficients of the unitary transform. In addition, we have omitted in Eq. (2b) the terms 
<4- c Lt an d d a Ck a , because they do not conserve the number of particles and could lead only 
to insignificant small corrections in comparison with a contribution of the term (2b). The 
bare vacuum of the model is defined now by Ck CT |0) = d a \0) = 0. To obtain the ground state 
of the BCS Hamiltonian, one needs thus to fill all states with k < kp. 

Making use of the spherical harmonic representation for band electron operators [0, we 
arrive at the following one-dimensional form of the model Hamiltonian: 
dk 



H = £ f°° 7T- {u(k)cl(k)c a (k) + v \dlc a (k) + ci(k)d a ] } + e d £ S a d a + Ud^d^ 

+ | J — OO 2:71 ^ L J ^ + | 

(3a) 



°-=T,l " ~°° <r=U 



where 



|A| 2 + P ,k<0 
u{k) = { v l_ (3b) 

\A\ 2 + k 2 ,k>0 



and the function e(k) is linearized around the Fermi level, e(k) ~ Vp(k — kp). We also have 
set Vf equals to 1 and count the electron momentum k from kp and the energy from ep. 
In the normal state, Eqs. (3) are reduced to the integrable version of the Anderson model. 
Hereafter, we restrict ourselves to very large values of U, so that the double occupancy of 
the ci-level is unlikely and can be excluded. 

The model Hamiltonian (3) describes the behavior of a transition metal impurity in a 
superconductor. In the case of rare-earth impurities, one needs to combine the BCS model 
and the degenerate Anderson model After analogous manipulations, we then obtain 

— \uj(k)ci(k)c a (k)+v cl(k)X 0a + X a0 c a (k) } + e f J2X aa . (4) 
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Here, the Fermi operator c* a {k) creates an electron with the total angular momentum projec- 
tion a and the momentum modulus k. An impurity is described by the Hubbard operators 
X a b with algebra X a bX c d = 5 c bX ac i, where the index a = 0, a enumerates both the nonmag- 
netic state (a = 0) and the degenerate magnetic states (a = 1, . . . , n) of the impurity with 
the /-level energy 6f. The Coulomb repulsion on the impurity orbital is assumed to be very 
large, such the multiple occupancy is excluded. In the normal state, Eq. (4) is reduced to the 
n-fold degenerate Anderson model, which has been diagonalized by Schlottmann || within 
the LDA. In the particular case n = 2, the model (4) is equivalent to the nondegenerate 
model (3) with infinitely large Coulomb repulsion U. 

To diagonalize the model Hamiltonian, it is convenient to rewrite Eq. (4) in terms 
of electron operators on the "energy scale", c a (u) = (duj(k)/dk)~ 1 ^ 2 c a (k), with algebra 
{c a (uj), c\,{oj')} = 2iT5 aa >5(u - u'), 

} + e f J2X aaj (5) 

a 

where V(u) = v(dcu(k)/dk)~ 1 ^ 2 is energy dependent because of the nonlinear electron dis- 
persion (3b). The integration contour C consists of two semi-infinite intervals, C = C_©C + , 
where C_ = (— oo, —A) and C + = (A, oo) correspond to the lower (k < 0) and upper {k > 0) 
branches of the electron dispersion. 

We look now for one-electron eigenstates of the Hamiltonian (5) in the form 

l*i) = EA 

a 

where A a is arbitrary. The vacuum state contains no band electrons and the impurity is in 
the nonmagnetic state, c a (e)|0) = X aa \0) = 0. The Schrodinger equation (H — = 

then reads 

{-id T - u)<j>{r\u) + g(w)5(T) = 0, (6a) 

/oo 
drT(T)<P(T\u) = 0, (6b) 
-oo 

where 0(r) is the Fourier image of the auxiliary wave function <p(e), while the effective 
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cl(u)c a (uj) + V(u) cl(uj)X 0a + X a0 c a (uj) 



X a o + / c |^)0(e)ci 



|0>, 



particle-impurity coupling T(r) contains the information about the electron dispersion, 

^ = /II^< r(r) = / o £v (e)e -. 

For w 6 C, one gets g(u) = [r(u) + and 

,/ I \ r M - (V2)sgn(r) jUjjt m 



The "rapidity" r(u;) is defined by 



u — e 

r = 



(8) 



where \^ 2 (o;)/2 is the imaginary part of the self-energy 

em-/* y2(t) . 

v ; Jc 2tt e - w - iO 

The real part of the self-energy, = ReS(u;), equals to zero if cj lies outside the gap, 

| a; | > A. If <jj lies inside the gap, uj G G = (—A, A), the imaginary part of the self-energy 
vanishes, and Eqs. (6) admit also the discrete mode with the eigenenergy ua, which is 
found root of equation e/ — uj = E'(o;), where E'(u) = v 2 u/2VA 2 -u 2 and u e G. The 
discrete mode describes the electron-impurity bound state, which is a complete analog of the 



discrete mode predicted earlier by John and Wang [lOj for a Bose system. If e/ lies far from 
the gap, \ef\ ^> A, the bound state energy Ud is very close to the conduction band states, 
but becomes well separated if e/ lies inside and around the gap. Note also that for many 
superconducting magnetic alloys v 2 ^> A, therefore Ud lies close to zero (the Fermi energy 
of the host metal) for arbitrary position of the impurity level energy e/. It should be also 
emphasized that the discrete mode is found here as a solution of the one-particle problem 
(6) rather than a multiparticle bound state discussed earlier by Shiba |[L1|| . Therefore strong 
electron-electron correlations could result in a strong renormalization of the bare discrete 
mode energy ujd- 

The auxiliary wave function (7) is discontinuous, but the electron wave function 
if) a (x\u) = (c a (x)|\E'i) = ip(x\u))A a (where the operator c a (x) is defined as c a (x) = 



J^° 00 (dk/2'ir)c a (k) exp (ikx) [[|]) is continuous and results from the integral dressing of the 
auxiliary function, 

/oo 
dru(x; t)4>(t\u), 
-oo 

with the dressing function 

In the LDA, the dressing function is nothing but the delta-function, u(x; r) = v8(x — r), and 
hence the auxiliary wave function (J)(t\uj) is the same as the electron wave function ip(x\u). 

The idea of auxiliary functions (or "auxiliary particles") can easily be extended to the 
multielectron case. For instance, the two-electron wave functions in the energy space are 
represented as ^ aia2 (e 1 ,e 2 ) = V(e 1 )V(e 2 )^ aia2 {ex, e 2 ) and J QlQ2 (e) = V(e)G aia2 (e), where 
the latter describes the state, in which one of the electrons is localized on the impurity. 
In the auxiliary r-space, the Schrodinger equations for the Fourier images of the auxiliary 
functions Q ai a 2 (. T i, r 2 ) and G aia2 (r) are then similar to those in LDA but with the nonlocal 
particle-impurity coupling T(r). It is remarkable that, despite the nonlocal coupling, the 
two-particle scattering matrix is still found in the well-known form: 

S= r( " l) - r( " 2 \ +2P (9) 

where P = 5 ai a > 5 a2a i is the permutation operator. The S'-matrix satisfies obviously the 
Yang-Baxter equations. Hence the multiparticle scattering is factorized into two-particle 
ones and the auxiliary multiparticle wave functions have the ordinary Bethe ansatz structure. 
But due to the nonlinear dispersion, the multielectron wave functions 

roo n 

(xx,...,X N )= $ ai _ a2 (Tx,...,TN)'Y[u(x j ;Tj)dTj (10) 

J-oo j=1 

are continuous functions of the coordinates Xj. The factorization of multielectron scattering 
and the Bethe ansatz construction for multielectron wave functions are thus hidden and 
manifested only in the limit of large interelectron separations. 
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Imposing the periodic boundary conditions on the iV-electron wave function (10) on the 
interval of size L, we arrive at the following Bethe ansatz equations (BAE): 

Mi 

exp (i%L)e x (Af ) = n ei(Af - A«) (11a) 

ai=l 

M l+V Mi 

n n ^-AK^-n^-A?) (n b ) 

^=±1 a; + „=l b=l 
AT 

E = 5>,- (lie) 

where E is the eigenenergy, tUj = ui(kj) is the energy of a "charge" excitation of the system 
with the momentum kj, and e n (x) = (x — in/2) / (x + in/2) . If m a is the number of electrons 
with the "color" a, the numbers Mi are defined by Mi = Y^=i m a, Mo being the total number 
of electrons, Mq = N. It is clear that only charge excitations with rapidities = r(kj) 
contain the information about the nonlinear dispersion, while BAE for the color rapidities 
{A^}, I — 1, . . . , n — 1 coincide with the corresponding equations in the LDA. 

In the particular case n = 2, Eqs. (11) are reduced to the BAE of the nondegenerate 
model (3) with infinitely large U, 

exp (ik.L) rj - 1,2 = n r i- X «- 1 / 2 (12 a ) 

nra-nra (12b) 

where the spin projection of the system is given by S z = N/2 — M. In the normal state, 
where uj = k and r(k) = (k — €f^)/v 2 , Eqs. (11) and (12) are the same as the well-known 
BAE of the Anderson model. 

In conclusion, it should be emphasized that the approach developed in the present paper 
for superconducting magnetic alloys is easily generalized to Fermi systems with nonlinear 
spectrum of particles around the Fermi level. For instance, the BAE derived here are valid 
also for gapless Fermi systems pi, where the function r(k) takes the form 



r{k) = "W-'z-Mq (13) 

v z dk 



and u(k) is the dispersion of particles of the system. 
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